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The Cluster-Merger Shock in IE 0657-56: Faster than a 
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Abstract. The merging galaxy cluster IE 0657-56, known as the "bullet clus- 
ter," is one of the hottest clusters known. The X-ray emitting plasma exhibits 
bow-shock-like temperature and density jumps. The segregation of this plasma 
from the peaks of the mass distribution determined by gravitational lensing has 
been interpreted as a direct proof of collisionless dark matter. If the high shock 
speed inferred from the shock jump conditions equals the relative speed of the 
merging CDM halos, however, this merger is predicted to be such a rare event in 
a ACDM universe that observing it presents a possible conflict with the ACDM 
model. 

We examined this question using high resolution, 2D simulations of gas 
dynamics in cluster collisions to analyze the relative motion of the clusters, the 
bow shock, and the contact discontinuity, and relate these to the X-ray data for 
the bullet cluster. We find that the velocity of the fluid shock need not equal 
the relative velocity of the CDM components. An illustrative simulation finds 
that the present relative velocity of the CDM halos is 16% lower than that of 
the shock. While this conclusion is sensitive to the detailed initial mass and 
gas density profiles of the colliding clusters, such a decrease of the inferred halo 
relative velocity would significantly increase the likelihood of finding IE 0657-56 
in a ACDM universe. 



1 . Int ro duct ion 

We use gas dynamical simulations of the nearly head-on collision of two unequal- 
mass clusters to model the "bullet cluster" IE 0657-56, to show t hat the high 
shock velocity of 4700 km s -1 infer red from X-ray observations (iMarkevitchl 
20061 : iMarkevitch fc Vikhlininl [2007) may exceed the relative velocity of the 



merging halos. Previous estimates of the extremely small probability of find- 
ing such a high-velocity cluster-cluster collision in a ACDM u niverse assumed 



that t he shock and the halo -collision velocities were the same (Havashi & White 



2006: iFarrar fc Rosenl 120071 ) . Our results, first described in iMilosavljevic et al. 
(2007), significantly increase the likelihood of observing such a merger event in 



ACDM. 
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2. Simulation 



We u sed the adaptive mesh-refinement (AMR) ASC FLASH code ([Frvxell et al 
2000) in 2D cylindrical coordinates to simulate a region 20 Mpc across, with 
1 kpc spatial resolution. Gravity is contributed by two rigid, spherical DM 
halo profiles, whose centers free-fall toward each other from rest at separation 
4.6 Mpc. The density associated with the gravitational force is an NFW pro- 
file p = /?o(r/r s ) _1 (l + r/r s )~ 2 for the subcluster, and a cored profile p = 
po(l + r/r s )~ 3 for the main cluster, to imitate a non-head-on merger. The 
masses (M^qq / 10 15 M Q ) and radii rsoo and r s (Mpc) are (1.25,1.5,0.5) and 
(0.25,0.5,0.16) for the main cluster and subcluster, respectively. 




Figure 1. Time evolution of temperature. Larger CDM halo is centered at 
cylindrical (r, z) — (0,0); smaller halo is at z — 0.29, —0.66 and —0.96 Mpc, 
respectively. Circles indicate the scale radii r s of the halos. 



3. Results 

After subhalo "bullet" passes through the center of primary halo, the subhalo 
gas trails the subhalo DM, led by a bow shock and contact discontinuity (cold 
front)., as seen in Figure U The opening angle and radius of curvature of the 
bow shock are sensitive to simulation details, but both are larger than expected 
for steady-state bow shocks of hard spheres with constant velocity in a uniform 
medium. The wings of the contact discontinuity are Kelvin-Helmholtz unstable. 

During pericenter passage, the shock and contact discontinuities are slower 
than the relative velocity of the two CDM halos, due to ram pressure force acting 
on the gas. Later, the halos climb out of the gravitational potential well and 
decelerate, but the shock and contact discontinuity do not decelerate (Fig. [2[a]). 
At observed separation D = 720 kpc, the velocity of the shock is 4800 km s _1 , 
consistent with the 4740 km s _1 inferred from IE 0657-56 data. The relative 
velocity of the simulated halos then is much less, 4050 km s _1 (16% smaller). 
The X-ray surface brightness and temperature profiles of the simulation viewed 
transverse to collision axis (Fig. |3j) agree with those of IE 0657-56. 
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Figure 2. (a) Velocity of the 
shock, "K;hock (solid line), and 
of the contact discontinuity, 
Contact (dashed line) relative 
to the preshock gas upstream, 
and the relative velocity of 
the CDM halos, V 12 (dash-dot 
line), all as functions of the 
time. The distance between 
two CDM halos, D, is used as 
the coordinate of time. (b) 
Density and temperature jump 
at the shock, and the shock 
Mach number, (c) Position of 
the shock (solid line) and the 
contact discontinuity (dashed 
line) relative to subhalo center. 
Figures [2] and [3] are taken from 
iMilosavlievic et alJ (|2007f ). 
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Figure 3. (Left:) Surface brightness map of the simulated cluster in the 
0.8 — 4 keV band. The circles indicate the scale radii r s of the CDM ha- 
los. (Middle:) Temperature profile along collision axis (solid) and emissivity- 
weighted temperature, averaged over the (transverse) line of sight (dotted). 
(Right:) Surface brightness profile along collision axis. 
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Figure 4. Thermal bremsstrahlung 
emission spectrum from simulated 
merging cluster: the contribution 
from cold gas below 15 keV (dotted) 
and hot gas above 15 keV (short- 
dash), total spectrum (solid), which is 
the sum of cold and hot components, 
and single-temperature fit to cold 
(long-dashed) and hot (mixed dash) 
components. The emission from the 
hot shock- heated gas (T > 15 keV) 
dominates the spectrum above 20 keV 
and is fitted well by a single tempera- 
ture profile with T = 20 keV. 



Figure H shows the thermal X-ray spectrum of the simulated merging clus- 
ter. The observed spectrum of IE 0657-56 has str onger emission above 30 keV, 



suggesting some additional non-thermal emission (|Petrosian et al.l 12006). 



4. Discussion and Conclusions 



Our simulations of merging galaxy clusters show that the halo collision velocity 
need not be the same as the intergalactic gas shock velocity. Applying our cor- 
rection to the speed of the subhalo, the probability that such a massive cluster in 
our CDM universe has such a high-speed subhalo increases by 3 orders of magni- 
tude relative to previous estimates, but is still small (~ 10~ 4 ). Springel & Farrar 
(j2007l ) independently reached a similar conclusion, although their relative halo 
velocity is much smaller, 2600 km s _1 . They assumed perfectly head-on collision 
between a cuspy main halo and a more concentrated subhalo than ours, which 
would make the subhalo deceler ate more before reaching the ob served separa- 
tion. More recent simulations by iMastropietro Sz Burkertl (|2007l ) show that the 
difference between the shock velocity and the halo relative velocity is smaller 
if the initial halo collision velocity is larger and if the collision is not perfectly 
head-on. 
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